Study of QCD critical point using canonical ensemble method 



Anyi Li^ Andrei Alexandru'' Xiangfei Meng and Keh-Fei Liu^, for the ;tQCD Collaboration 



The existence of the QCD critical point at non-zero baryon density is not only of great interest for 
experimental physics but also a challenge for the theory. We use lattice simulations based on the 
canonical ensemble method to explore the finite baryon density region and look for the critical 
point. We scan the phase diagram of QCD with three degenerate quark flavors using clover 
fermions with m„ ^ 700MeV on 6^' x 4 lattices. We measure the baryon chemical potential as 
we increase the density and we see the characteristic "S-shape" that signals the first order phase 
transition. We determine the phase boundary by Maxwell construction and report our preliminary 
results for the location of critical point. 



1. Canonical ensembles 

The search for the QCD critical point has attracted considerable theoretical and experimental 
attention recently. To simulate matter at high density, a method based on the canonical parti- 
tion function was proposed [IJ. While expensive - every update involves the evaluation of the 
fermionic determinant - finite density simulations based on this method proved feasible |2|. 

We construct canonical partition function as a Fourier transform of grand canonical partition 
function with respect to a t/(l) phase in the last time slice: 



where detitM^(U) is the projected determinant on the k quark sector. With the aid of winding 
number expansion method |3 1, this can be computed efficiently. A program was outlined to scan 
the QCD phase diagram to look for the critical point |4 5|. We refer readers to our previous 
papers that include details on the simulations of the canonical partition function |2 |. 

In this paper, we present results for Nf - 3 based on simulations on 6^' x4 lattices with clover 
fermions. We fix temperature and scan in the baryon number direction. By taking the difference 
of free energy after adding one baryon, baryon chemical potential can be measured as an observ- 
able in canonical ensemble. We plot the chemical potential as a function of baryon density. In 
a finite volume, due to the non-zero contribution from the surface tension, the first order phase 
transition will be reflected as an "S-shape" structure in this plot l6|. The phase boundaries of the 
coexistence region can be determined by "Maxwell construction". We observe a clear signal for 
the first order phase transition. With the phase boundary determined at simulation temperatures 
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Figure 1: "S-shape" structure in tlie baryon cliemical potential vs. baryon number plane as well as Maxwell construction 



by "Maxwell construction", we located the critical point at the intersection of the extrapolated 
phase boundary lines. 



2. Results 

We illustrate the "S-shape" structure in Fig.[T]at a fixed temperature T - O.SSTc. Once we 
determine the phase boundaries at a few temperatures, the critical point can be located at the point 
where phase boundaries of coexistence phase cross together Using the Maxwell construction we 
also determine the value of the critical chemical potential and we plot the phase diagram in the 
T - hb plane. 




Figure 2: Left panel: Boundaries of coexistence region. Right panel: Phase transition line in the T ~ plane. The 
critical point is located by an extrapolation 

The critical point is found at Te/Tc - 0.93(2) and H^IT^ - 3.25(7) for Nf - 3 and m„ » 
700 MeV. 
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